Instability of nuclear wobbling motion and tilted axis rotation 

Masayuki MatsuzakQ and Shin-Ichi Ohtsubo 

Department of Physics, Fukuoka University of Education, 
Munakata, Fukuoka 811-4192, Japan 
(Dated: February 5, 2008) 

"^j-" ■ Abstract 
O 

We study a possible correspondence between the softening of the wobbling mode and the "phase 

CN , 

transition" of the one-dimensionally rotating mean field to a three-dimensionally rotating one by 

a • 

comparing the properties of the wobbling mode obtained by the one-dimensional cranking model 
■ + random phase approximation with the total routhian surface obtained by the three-dimensional 

tilted-axis cranking model. The potential surface for the observed wobbling mode excited on the 

(N 

' triaxial superdeformed states in 163 Lu is also analyzed. 
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I. INTRODUCTION 



The concept of the phase transition of the mean field is useful for describing structure 
changes in the atomic nucleus although it is a quantum system composed of finite num- 
ber of Fermions. Typical example is that a spherical mean field becomes unstable as the 
quadrupole vibration excited on top of it softens with changing particle numbers, then an 
axially symmetric mean field substitutes. This can rotate about one of the axes perpen- 
dicular to the symmetry axis. Consecutively, the axial symmetric mean field can become 
unstable as the 7 vibration softens, then a triaxially deformed mean field substitutes. This 
can rotate about all the three principal axes. Usually, however, a rotation about one axis 
dominates because the rotation about the axis with the largest moment of inertia is ener- 
getically favorable. When some excitation energy is supplied, small rotations about other 
two axes become possible. Consequently this produces a kind of vibrational motion of the 
rotational axis, that is, the wobbling motion. 

The small amplitude wobbling motion at high spins was first discussed by Bohr and 
Mottelson in terms of a macroscopic rotor model. Then it was studied microscopically by 
Janssen and Mikhailov [2] and Marshalek [3] in terms of the random phase approximation 
(RPA). Since the small amplitude wobbling mode has the same quantum number, parity 
7r = + and signature a = 1, as the odd-spin member of the 7 vibrational band, Mikhailov 
and Janssen 4] anticipated that it would appear as a high-spin continuation of the odd-spin 
7 band. But it has not been clear in which nuclei, at what spins, and with what shapes it 
would appear. Using the RPA, Shimizu and Matsuyanagi |5| studied Er isotopes with small 
I7I, Matsuzaki (| and Shimizu and Matsuzaki 7] studied 182 0s with a rather large negative 
7 but their correspondence to the experimental data was not very clear. In 2001, Cdegard 
et al. jsj found an excited triaxial superdeformed band in 163 Lu and identified it firmly as 
a wobbling band by comparing the observed and theoretical interband E2 transition rates. 
These data were investigated in terms of a particle-rotor model by Hamamoto |S| and in 
terms of the RPA by Matsuzaki et al. In 2002, two-phonon wobbling excitations were 
lls0 „« by w e u ( .fl-l al excitation energies show some anharmonicity. 

The one-dimensionally rotating triaxial mean field may become unstable as the wobbling 
mode softens with changing some parameters. One of the present authors (MM) pointed out 
its theoretical possibility in the preceding Rapid Communication 10]. The possibility of this 
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phase transition was discussed in terms of the harmonic oscillator model by Cuypers |l2j and 
Heiss and Nazmitdinov [13j but their conclusions are controversial. A theoretical framework 
to describe three-dimensional rotations, possibly with large amplitude, was first devised by 
Kerman and Onishi 14] within a time-dependent variational formalism. Onishi Q| and 
Horibata and Onishi 16 1 applied it to 166 Er and 182 0s, respectively. See Ref. for exam- 
ple, for recent applications. The three-dimensional cranking model was first used by Frisk 



and Bengtsson 



18| . The word, "tilted (axis) cranking (TAC)" was, to our knowledge, first 



used by Frauendorf ^| and it was applied to a kind of two-dimensional rotation — the 



so-called shears band, observed for example in the A ~ 200 region |20j, |21 
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22J. Applications 
and references 



to multiquasiparticle high-i^ bands were also extensively done; see Ref. 
cited therein. When the rotation becomes fully three-dimensional, the new concept, chi- 
rality, emerges Q Q. The tilted mS caking was also applied to this Q. At finite 
temperature, the degree of freedom of spin orientation was studied macroscopically [27[ 
and microscopically |2£ LA relativistic formulation of the three-dimensional cranking was 
given by Kaneko et al. [29] as an extension of the one-dimensional one given by the Munich 
group 0| . Madokoro et al. 0] studied the shears band in 84 Rb starting from the meson 
exchange interaction although the pairing field was neglected. 

The purpose of the present paper is to elucidate the work in Ref. |l0j| by comparing two 
types of theoretical calculations, the one-dimensional cranking model + RPA and the three- 
dimensional (tilted axis) cranking model. The former gives the "mass parameters" for the 
motion of the angular momentum vector, that is, the moments of inertia, while the latter 
provides the surfaces on which the angular frequency vector moves around. 



II. MODEL 



We start from a one-body Hamiltonian in the rotating frame, 



h! = h — h cr , 
h = h mi -A T {P} + P T )-\ T N T , 

hmi = ^ + 2 M (^ 2 + ^ + w ^ 2 ) 
+^J-s + ^(l 2 - (l% osc ). 



(1) 
(2) 

(3) 



3 



{her is specified below.) In Eq.(P|), r = 1 and 2 stand for neutron and proton, respectively, 
and chemical potentials A T are determined so as to give correct average particle numbers 
(N T ). The oscillator frequencies in Eq.(jSJ) are related to the quadrupole deformation pa- 
rameters €i and 7 in the usual way. (We adopt the so-called Lund convention.) They are 
treated as parameters as well as pairing gaps A T . The orbital angular momentum 1 in Eq.(jSJ) 
is defined in the singly-stretched coordinates x' k = y^^x^, with k = 1 — 3 denoting x — z, 
and the corresponding momenta. 



A. One-dimensional cranking model + random phase approximation 

Equations (|IJ) - © with 

h cr = hhJ Iot J x (4) 

generate the system rotating one-dimensionally. Then, since h' conserves parity tt and 
signature a, nuclear states can be labeled by them. Nuclear states with quasiparticle (QP) 
excitations are obtained by exchanging the QP energy and wave functions such as 

(-e^V^U^-^.Up.Vp), (5) 

where p, denotes the signature partner of \i. We perform the RPA to the residual pairing 
plus doubly-stretched quadrupole-quadrupole (Q" ■ Q") interaction between QPs. Since we 
are interested in the wobbling motion that has a definite quantum number, a = 1, only two 
components out of five of the Q" ■ Q" interaction are relevant. They are given by 

rr(-) _ 1 ,,(-W'(-) W'H (R\ 

H int ~ -g K K Qk Qk > K b ) 

if =1,2 

where the doubly-stretched quadrupole operators are defined by 

Q'k = Qk{xic -> x'L = —x k ), (7) 
and those with good signature are 

Q ( k ] = /o / = f (Qk ± Q~k) ■ (8) 
a/2(1 + 6 KQ ) 

The residual pairing interaction does not contribute because P T is an operator with a = 0. 
The equation of motion, 

'h' + H^\xl\_ =?iu n Xl (9) 



J RPA 
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for the eigenmode 



(a=±l/2) 



pL<U 

leads to a pair of coupled equations for the transition amplitudes 
Then, by assuming 7 7^ 0, this can be cast H into the form 



(10) 



(11) 



^rot) 



2 2 
_ ^rot 



Kcflf) 



(eff), 



0. 



(12) 



which is independent of k^s. This expression proves that the spurious mode (u n = u> I0t ; 
not a real intrinsic excitation but a rotation as a whole) given by the first factor and all 
normal modes given by the second are decoupled from each other. Here J x = {J x )/uj rot as 
usual and the detailed expressions of Jy? 

(oj n ) are given in Refs. |3j, |y,|2{- Among normal 

modes, one obtains 



^wob — ^rot 



\ 



iyJx~ Jy ^(^wob)j {jx" Jz ^ 



^wob, 



(13) 



Jy \il) wo t)Jz ^(^wob) 

by putting u n = co> wo b- Note that this gives a real excitation only when the argument of 
the square root is positive and it is non-trivial whether a collective solution appears or not. 

n 

Evidently this coincides with the form derived by Bohr and Mottelson in a rotor model [1( 



and known in classical mechanics 



32 ]. 



B. Three-dimensional (tilted axis) cranking model 

In this model the one-body Hamiltonian is given by Eqs.Q - (JHJ) with 

h cr = • J, (14) 
f2 = dj rot (cos^, sin# cos <p, sin 6 sin ip). (15) 

Pairing correlation is taken into the system by a simple BCS approximation with fixed gaps 
as in the case of the one-dimensional cranking. The expectation value (J) calculated at each 
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(urot, 0, ip) has three non-zero components in general; the stationary state that minimizes 
the total routhian is obtained by requiring (J) || f2 (see Ref. for details). Obtained tilted 
solutions do not possess the signature symmetry and therefore describe AJ = 1 rotational 
bands. In the present work, given a set of mean-field parameters, N T , €2, 7, and A r , a 
configuration is specified at 9 = 0° (principal axis cranking about the x axis). Then by 
changing 6 and tp step by step, the most overlapped state is chased. This procedure gives 
an energy (total routhian) surface for the angular frequency vector. Surfaces for QP excited 
configurations can also be calculated by adopting a procedure similar to Eq.©. 



III. RESULT AND DISCUSSION 

For this first comparative calculation, we choose the [(^9/2, /7/2) 2 ( 7r ^n/2) 2 ]i6+ f° ur quasi- 
particle configuration in 146 Gd among this mass region in which many oblate isomers have 
been observed. This state is described by e 2 = 0.19, 7 = 60°, A n = 0.8 MeV, A p = 0.6 MeV, 
and tijjJrot = 0.25 MeV. Calculations are performed in the model space of three major shells; 
-Wosc = 4-6 for neutrons and 3-5 for protons. The strengths of the 1 • s and l 2 potentials 
are taken from Ref. 

In the present study we concentrate on the changes in the system with 7. Figure E^a) 
reports the excitation energy hw wo b in the rotating frame. That in the laboratory frame in 
the case of 7 = 60° is given by ftw wo b + fo^rot = 0.198 MeV + 0.25 MeV. The excitation 
energy decreases steeply as 7 decreases. In order to see its implication, we show in Fig. [Tfb) 
the wobbling angles, 

^IJ^Ko^P + IJi^Kob)! 2 

6 wob = tan" 1 ^ — , (16) 



V?wob = tan 1 



Jz (^wob) 



(17) 



Jy \^wob) 

with (PA) denoting the principal axis. # wo b clearly proves that the softening of the excitation 
energy is accompanied by a growth of the amplitude of the motion. p wo b indicates that the 
fluctuation to the y direction grows. Corresponding to this, the three moments of inertia 
behave as in Fig.^c). Qualitatively, this behavior can be understood as an irrotational-like 
moments of inertia, 

Jj? ocsin 2 ( 7 + ^), (18) 



where k = 1 - 3 denoting the x — z components, superimposed by the contribution from the 
alignment, AJ~ X . Alternatively, it can also be viewed as that, at large 7, multiple alignments 
lead to a rigid-body-like inertia, 



Ji ig oc 1 1 




J^/3cos (7 + ^irk) 



(19) 



with (3 being a deformation parameter defined by the mass distribution. 
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FIG. 1: Triaxiality dependence of (a) excitation energy of the wobbling motion, (b) wobbling 
angles, and (c) three moments of inertia associated with it in 146 Gd, calculated at hu TO t = 0.25 
MeV with e 2 = 0.19, A n = 0.8 MeV, and A p = 0.6 MeV. 



Now we proceed to three-dimensional calculations; we calculate energy surfaces as func- 
tions of the tilting angle (9, ip) of f2. Here we note that the (9, tp) plane is represented as a 
rectangle although ip is meaningful for 9^0°. Figure Efa) shows the 7 = 60° (symmetric 
about the x axis) case. Until down to 7 ~ 40°, energy surfaces are qualitatively similar aside 
from becoming shallow gradually. But a further decrease of 7 leads to an instability of the 
motion to the 9 direction with ip = 0°, that is, the direction of the y axis. Together with the 
property that the surface is stable with respect to the direction of the z axis, the situation 
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corresponds excellently to Fig. [T] The behavior of (p wo b in Fig. [TJb) can be interpreted as 
that, when the system can fluctuate to the direction of the y axis without any energy cost, 
it does not fluctuate to the z axis. 






FIG. 2: Energy surfaces of the [(W19/2, /7/2) 2 ( 7r ' i ii/2) 2 ]i6+ configuration in 146 Gd as functions of 
the tilting angle (9,<p) calculated with the same parameters as Fig. E ( a ) 7 = 60°, (b) 7 = 40°, 
(c) 7 = 30°, (d) 7 = 20°, and (e) 7 = 0°. The interval of contours is 50 keV. Discontinuities in the 
surfaces are due to quasiparticle crossings. 

To look at the energy surface more closely we gather their cross sections at tp = 0° (the 
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x-y plane) in Fig. EJ This figure clearly shows that a tilted axis minimum appears at around 
7 = 30° although it is shallow. The correspondence to Fig. in which the instability occurs 
at 7 = 32° is excellent. Note that the reason why the wobbling angle seen from Fig. |H]is 
larger than 8 wo ^ in Fig. ^b) is that this is drawn for f2. 
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FIG. 3: Cross sections at (p = 0° of the energy surfaces of 146 Gd. 



IV. POTENTIAL SURFACE FOR THE WOBBLING MODE IN 163 LU 

The analyses above are purely theoretical. Then, is there any experimental signature 
of the softening of the wobbling motion? We think the answer is yes. Figure H] shows 
the experimental excitation energies (in the rotating frame) of the TSD3 (two-phonon 
wobbling) and the TSD2 (one-phonon wobbling) relative to the TSD1 (yrast 1QP TSD) 
in 163 Lu, where TSD is the abbreviation for triaxial superdeformation. Ai?2_ p h on on < 2 x 
A_E(_ phonon indicates a signature of softening of the energy surface. We obtained hw wo b = 
0.185 MeV, flwnh = 14.2°, and y? wo b = 7.6° for the one-phonon wobbling state in the RPA 
(see also Refs.jiol for the RPA calculation). The small value of y? wo b looks to indicate a 
softening to the y direction. Calculated energy surface is shown in Fig. |S| Calculations were 
done in the model space of five major shells, A osc = 3-7 for neutrons and 2-6 for protons, 
with e 2 = 0.43, 7 = 20°, A n = A p = 0.3 MeV, and frw Iot = 0.5 MeV where the calculated 
^wob approaches the experimental one. This figure shows again the surface softens to the 



direction o 
In Refs. 



the y axis. 

lOj, |34j, it was shown that the alignment of the 7ri 13 / 2 quasiparticle was essential 



for the appearance of the wobbling motion. In order to see this fact from the viewpoint of 
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FIG. 4: Experimental excitation energies of the two- and one-phonon wobbling states relative to 



the yrast triaxial superdeformed states in 163 Lu. Data are taken from Ref. 








FIG. 5: Energy surface of the triaxial superdeformed one-quasiparticle configuration in 163 Lu as 
a function of the tilting angle (0,(f) calculated at Huj ro t = 0.5 MeV with e 2 = 0.43, 7 = 20°, and 
A n = A p = 0.3 MeV. The interval of contours is 100 keV. 

the potential surface, we calculated the 0QP (non-yrast TSD at high spins) configuration 
in 162 Yb. Figure EH clearly shows that a tilted axis minimum in the x-y plane is realized 
when the wobb ling motion does not occur due to the lack of vri 13 / 2 QPs that make J x larger 
than jiol 34| . This result proves that the low-f2 high-j orbital favors the principal axis 

rotation on which the wobbling motion occurs. 

For deeper understanding of the two-phonon states, application of more sophisticated 
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FIG. 6: The same as Fig. but for the zero-quasiparticle configuration in 162 Yb. 

many body theories such as the selfconsistent collective coordinate (SCC) method is 
desirable. 

V. SUMMARY 

To summarize, we have proved that a tilted axis rotation emerges when the wobbling mode 
becomes unstable as the triaxiality parameter changes in an oblate configuration in 146 Gd. Its 
instability is caused by a growth of the fluctuation of the motion of the angular momentum or 
frequency vector to the direction of the y axis. Having performed this theoretical calculation, 
we have argued that the signature of the softening of the wobbling motion can be seen in 
the observed spectra of the triaxial superdeformation in 163 Lu and shown that a tilted axis 
minimum would appear if it were not for the 7ri 13 /2 quasiparticle. 
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